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Abstract 

The mathematical consequences of a few simple scaling assumptions about the effects of compress- 
ibility are explored using a simple singular perturbation idea and the methods of statistical fluid 
mechanics. Representations for the pressure-dilatation and dilatational dissipation covariances 
appearing in single-point moment closures for compressible turbulent are obtained. The results 
obtained, in as much as they come from the same underlying diagnostic relationship, represent 
a unified development for both the compressible covariances. While the results are expressed in 
the context of a second-order statistical closure they provide some interesting and very clear phys- 
ical metaphors for the effects of compressibility that have not been seen using more traditional 
linear stability methods. In the limit of homogeneous turbulence with quasi-normal large scales 
the expressions derived are - in the low turbulent Mach number limit - asymptotically exact. The 
expressions obtained are functions of the rate of change of the turbulence energy, its correlation 
length scale, and the relative time scale of the cascade rate. With the appearance of the length 
scale the dilatational covariances are found to scale with the Mach numbers based on the mean 
strain and rotation rates. The expressions for the dilatational covariances contain constants which 
have a precise and definite physical significance; they are related to various integrals of the longi- 
tudinal velocity correlation. The pressure-dilatation covariance is found to be a non-equilibrium 
phenomena related to the time rate of change of the internal energy and the kinetic energy of the 
turbulence. Also of interest is the fact that the representation for the dilatational dissipation in a 
turbulence, with or without shear, features a dependence on the Reynolds number. This article is 
a documentation of an analytical investigation of the implications of a pseudo-sound theory for the 
effects of compressibility. The novelty of the analysis is in the very few phenomenological assump- 
tions required to produce the results. Subsequent work will assess the consequences of this analysis 
in the context of compressible turbulence models for engineering calculations. 

*This research was supported by the National Aeronautics and Space Administration under NASA Contract No. 
NAS1-19480 while the author was in residence at the Institute for Computer Applications in Science and Engineering 
(ICASE), NASA Langley Research Center, Hampton, VA 23681-0001. 
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1 Introduction 

In the moment equations for compressible turbulence several new quantities, not seen in the in- 
compressible form of the equations, appear. There are new terms in the equations that reflect 
the fact that the fluctuating dilatation, d = in a compressible turbulence is not zero. These 
effects have been divided into two categories, Lele(1994). There are contributions to the fluctuating 
dilatation by the the fluctuating pressure; there are also contributions to the fluctuating dilatation 
by fluctuations in composition or temperature which occur in situation in which of heat and mass 
transfer are relevant. These two sources of the fluctuating dilatation have been distinguished using 
the phrases compressibility effects and variable inertia effects, Lele (1994). There are also additional 
effects, in inhomogeneous or nonequilibrium flows, associated with the portion of density fluctua- 
tions due to gradients in the mean density; these effects might also be called variable inertia effects 
giving rise to the differences between Favre and Reynolds averaged variables. There are also effects 
associated with the transport coefficient variations due to the fluctuations in temperature. In this 
article effects due the occurrence of a nonzero fluctuating dilatation, the so-called compressibility 
effects, are treated. 

In the context of single-point moment closure methods, compressibility effects due to the fluctuating 
divergence appear in two new terms in the kinetic energy equation of a turbulent field: the pressure- 
dilatation, < pd >, and the variance of the dilatation, < dd >, which is related to what has come to 
be called the compressible dissipation, € c = < dd >. The turbulent energy equation is written 

as 

< p > j^ k = P k - < p>e s + < pd > - < p > € c + T k . (1) 

P k represents the production and T k represents the transport terms and any other terms (that 
are not directly germane to the present analysis). There are additional terms representing the 
contraction of the mass flux vector on the mean flow acceleration. T k will also be used to represent 
all such terms. The < pd > and < dd > appearing in the k equation were recognized by Zeman 
(1990, 1991) and Sarkar et al. (1991) in earlier studies of compressible turbulence closures. They 
are the subject of this work. 

The dilatational covariances also appear in the internal energy equation, here written in terms of 
the mean temperature with a constant c v : 

< p > c v j^ T = Pj - < pd> + < p> € s + < p > e c + Tj 

Where Tj is the transport of the mean temperature including such effects as heat flux and 
the turbulent or pressure transport. The production, for a specific class of flows, is given by 
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Pj = -PD + 2 < fi > (5 2 + W 2 ) where the strain and rotation tensors are defined to be 
symbolically equivalent to the incompressible case, ie. traceless: =_ 2 [Viij + ^jn 3 Dfiij], 

Wij — — Ujyi]. Note that Sjj = 0 since D — Uj,j . There are additional terms depend- 

ing on heat or species transfer and fluctuations in fluid properties: they are not germane to the 
present analysis. Note that the dilatational covariances appear with opposite signs in the kinetic 
and mean internal energy equations. The dilatational covariances represent an irreversible, £ c , and 
reversible transfers, < pd >, of energy between the mean internal energy field and the fluctuating 
kinetic energy field. 

The pressure-dilatation covariance, < pd >, and the dilatational variances, < dd >, - which is 
related to the dilatational dissipation - have been the subject of several studies stressing both the 
fundamental issues in understanding the physics as well as obtaining models suitable for use in 
engineering calculations. In addition to the aforementioned works of Zeman and Sarkar, additional 
insight into these terms can be found in the studies of Durbin and Zeman(1992), Zeman and 
Coleman (1991), Zeman (1993), Erlebacher et ai (1990), Sarkar (1992), Sarkar et ai (1991a, 
1991b), Blaisdell et ai (1991), Blaisdell and Sarkar (1993), Lee (1992). 

The present approach differs from the approaches of both Zeman and Sarkar. A low turbulent Mach 
number expansion of the equation of state, the Navier Stokes, the continuity and wave equations is 
conducted. The problem is recognized as a singular perturbation in as much as there are two relevant 
length scales: an inner scale, £, associated with the turbulence field, and an outer scale A r C/Mt 
associated with a propagating “acoustic” radiation field surrounding the vortical motion producing 
the radiation field. The perturbation development produces an algebraic constitutive equation for 
the fluctuating dilatation: the continuity equation, rather than being a prognostic equation for the 
density, becomes a diagnostic equation for the fluctuating dilatation. Taking the relevant moments 
of the expression produces constitutive relations for < pd > and < dd >. Assuming homogeneity 
and quasi-normalitv expressions without any undefined constants are obtained for < pd > and 
< dd >. 

Retaining only the lowest order isotropic contribution produces simple expressions in terms of an 
incompressible turbulence for the unknown covariances are found. The expressions are then given 
in terms of quantities carried in a single-point closure. 

The analysis can be contrasted to more traditional approaches using linear stability theory. Lele 
(1994) provides a resume of several of these works. Despite the inherent limitations in the linear 
stability analysis there has been some useful light shed on the dynamical aspects of the effects of 
compressibility on the flow. Here a statistical approach to the problem is taken. Such an approach 
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accounts for the nonlinearity of the phenomena. A statistical approach approach, however, does 
average over the many interesting dynamical features of the fluctuation flow_but delivers some very 
interesting insights indicating the nature of the cumulative effects of the fluctuations. Of particular 
note is the interpretation of the effects of compressibility as an added mechanism for the transfer 
of energy between the turbulence field and the mean internal energy field. The theory predicts the 
mechanism and the rate of this intermodal energy transfer. 

It is found, in analogy with the pressure-strain covariances in the Reynolds stress equations for an 
incompressible flow, Ristorcelli et al (1994), that the representations for the dilatational covari- 
ances have a rapid and a slow component. However, unlike the rapid-pressure in the incompressible 
problem, the rapid component depends on the spatial area of the rapid-pressure correlation. This 
behavior results in a dependence on a “gradient” Mach number: a Mach number formed using 
the mean velocity gradient and a length scale of the turbulence. This quantity appears to be an 
important parameter in distinguishing the effects of compressibility in a mixing layer from those in 
a boundary layer, Sarkar (1994). 

Several other interesting results have been found. These are now highlighted: 

1) The pressure dilatation < pd > can be either positive or negative depending on the rate of 
change of the kinetic energy, the mean temperature, the length scale of the turbulence and the mean 
velocity gradients. For a near equilibrium flow, as long as the production exceeds the dissipation 
by an amount that scales with the square of the turbulent Mach number, M t 2 = — 3^, and the rate 
of increase of the internal energy field, the pressure-dilatation will be negative transferring energy 
from kinetic to internal modes. 

2) The representation for < pd > can be shown to behave as an added mass term in the the k 
equation: inertia is added to the turbulence by the capacitance of the fluctuating pressure field. 
The fluctuating pressure field, or equivalently the mean internal energy, acts as a capacitor storing 
energy fed into the turbulence by the production and then transferred to the internal energy by 
< pd >. 

3) The pressure dilatation, < pd >, scales with Mf with order M 4 and higher corrections. The 
dilatational dissipation scales as A/ 4 with a ReJ 1 dependence; R t = is the turbulent Reynolds 
number. For high Reynolds number flows, much higher than those seen in DNS, the dilatational 
dissipation is found to be small; the primary effects of compressibility are due to the pressure- 
dilatation. The effects of compressibility, occurring through the agency of the pressure-dilatation 
covariance, are found to be important in nonstationarv flows. 
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4) The rapid portions of both < pd > and < dd > scale with the relative time scales of the cascade, 
Sk/e and Wk/e , (the mean strain and rotation) to the second and fourth powers. With the 
concomitant appearance of the turbulent Mach number M t , the expressions are found to scale with 
the mean deformation and mean rotation Mach numbers Sl/c and Wife. It is these Mach numbers 
that distinguish the mixing layer from the equilibrium boundary layer and thus the representation 
can discriminate between these two types of flows. The appearance of the two different gradient 
Mach numbers indicates a dependence of the dilatational covariances on the relative amounts of 
strain versus rotation; a fact not yet looked for in experimental or numerical results. 

5) Unlike the incompressible dissipation, in which the viscosity sets the small scales but otherwise 
has very little effect on the cascade rate, the dilatation dissipation is, in this linear theory, dependent 
on the viscosity. For a fixed M t number, as R t — ► oo, the dilatational dissipation vanishes. The 
compressible dilatation doesn’t appear to be able to be understood as a spectral cascade rate set 
by the large scales of the flow. 

6) It is seen that the density fluctuations are related to the incompressible pressure fluctuations, 
p = ip which, upon rescaling, can be understood as, p = Mfp. This indicates that compressible 
numerical simulations starting from “incompressible” initial conditions are more consistently ini- 
tialized with nonzero initial density and temperature fluctuations. One can speculate that initial 
conditions inconsistent with the variances associated with the incompressible pressure may create 
a wave field that may delay the decay of the transients or that may not even decay during the 
course of a DNS. These transients are analogous to the transients associated with the free motion 
of a second order system, y" + y' + y = 0 relaxing from some initial condition. This situa- 
tion appears to correspond to the analysis followed by Erlebacher et al. (1990) and Sarkar et al. 
(1991b). The dependence on the initial conditions in compressible isotropic simulations has been 
thoughtfully noted in Blaisdell et al. (1993). The crucial point, as is made clear by the analysis, is 
that the initialization of any calculation with zero fluctuating temperature, density or dilatation is 
inconsistent with a finite non-zero turbulent Mach number. 

The present treatment for the effects of compressibility can be thought of as analogous to the forced 
system, y" + y' + y = /(cut); the forcing coming from the vortical motions of a turbulence with 
non-zero Mach number in which the effects of transients from initial conditions has faded. The 
analogy can be made exact, however, doing so is not relevant to the present subject. 

The present article is organized in the following fashion: governing equations, analysis, discussion of 
physic of the results of the analysis, discussion of limitations and assumptions. The first two sections 
are fundamental in laying the ground work for the representations for the covariances with the 
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fluctuating dilatation: first a simple heuristic picture of the physics is presented after which a system 
of equations consistent with the physics presumed is derived. In the subsequent section and its 
four subsections, the assumptions of homogeneity and isotropy are exploited in obtaining analytical 
expressions for the desired covariances. The methods of statistical fluid mechanics, following the 
inceptional works of von Karman and Howarth (1938), Batchelor (1951), and Proudman (1952), 
are relied on extensively. As a byproduct of the section on the rapid pressure-dilatation correlation 
an expression for the pressure variance in an arbitrary three-dimensional mean flow is derived. 

Later sections discuss the physical implications of the representations derived. Qualitative com- 
parison is made to the physics that is known for several simple compressible turbulent flows. It is 
shown that the models show no effects of compressibility in the equilibrium adiabatic wall layer as 
is known to be the case. The document finishes with a summary of the limitations and assump- 
tions built into the theory which suggest future work as well as the class of flows for which the 
representations are expected to be useful. 

This article is meant to be primarily analytical. A simple perturbation analysis and the methods 
of statistical fluid mechanics are used to investigate the implications of a few simple and reasonable 
assumptions. The results are a mathematical consequence of the initial assumptions. The article 
is intended to be a documentation of this procedure and its implications. The objective is to 
providing metaphors and nondimensional numbers with which to understand and further explore 
diverse issues in compressible turbulence. Testing, verifying, exploring and evolving the present 
analytical results into a working turbulence model suitable for engineering calculations is the subject 
of a sequel work now in progress. These more quantitative issues are addressed in several works 
planned and in progress, Ristorcelii (1995), Ristorcelli et al. (1995). 

2.1 A physical background for the mathematics 

Before presenting the mathematical development leading up to the analytical expressions for pres- 
sure dilatation and the dilatational covariances a physical picture underlying and suggesting the 
mathematical development is described. A more formal and mathematical presentation is given in 
due course. 

It is useful to keep in mind the one essential and central bit of physics that forms the lynchpin of the 
theory and makes the present method and results possible: in the near field of an acoustic source, 
whose size is small with respect to the wavelength of its emission, the fluid behaves as if it were 
incompressible. This observation appears to have been first made by Landau and Lipschitz (1958) 
and is a cornerstone in the method of matched asymptotic expansions in the field of acoustics. 
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Three basic ideas form the foundation of the pseudo-sound theory for the dilatational covariances. 

The first is recognizing the problem as a singular perturbation problem - it has two different length 
scales. As it is the turbulence that creates the pressure and density fluctuations in the medium 
the frequencies of the compressible disturbances are the same as the frequencies of the turbulence. 
c/A ~ iijt. The two length scales are: a correlation length scale associated with the fluctuations 
of the turbulence, l, and a length scale A ~ i/M t associated with the propagation of pressure 
and density fluctuations associated with the turbulent fluctuations. Here M t = (2k/3) 1 ^ 2 /c is 
the turbulent Mach number where k = 1/2 < ujitj > and c ^ = ~fPoo/p<x, is the sound speed. 
The turbulent Mach number is to be used as the small parameter in expansions of the compressible 
equations to obtain representations for the effects of compressibility as manifested in the covariances 
with the fluctuating dilatation. Underlying the low M t assumption which leads to the two disparate 
scales is what is called, in the sound generation problem of aeroacoustics, the compact-source 
assumption. This is equivalent to the idea that the flow structure covers a distance small with 
respect to the length scale of the compressible radiation it emits. Closely related to these two 
length scales are two time scales: one associated with the convective modes of the flow, sav u/£, 
and the sound crossing time - the time it takes for a information to cross a typical scale of the 
turbulence, c/t. Note that the conventional definition of the Mach number is used: it is the ratio 
of a characteristic fluctuating velocity to the (mean) sound speed. This is in concordance with the 
conventions of the acoustics literature from which some of our ideas are drawn. 

The second idea concerns the pressure. In the sound generation problem, two pressures, an “acous- 
tic” pressure which propagates and a pseudo-pressure associated with the convective motions of 
the fluid, are sometimes distinguished. The term pseudo-pressure was first coined by Blokhintsev 
(1956) as quoted in Ribner (1962). The term propagating pressure will be used for the term “acous- 
tic” pressure so as not to imply that the problem is linear as the propagation of sound is assumed to 
be in the small disturbance limit. Without attempting to be precise, the pressure fluctuations in a 
fluid satisfy, to the degree suitable to the present heuristic discussion, the following wave equation, 
Light hill (1952) 


*"oo Pitt Pijj ~ (pWjUj ),ij (2) 

where p represents the deviations of the fluid pressure from its value at the static reference state. 
Care must be taken in assessing which solutions to this equation are relevant to compressible tur- 
bulence modeling. Solutions to this equation are comprised of the homogeneous solution, which 
obeys the sourceless wave equation and the wave equation with source due to the turbulent fluc- 
tuations. The sourceless wave equation, essentially the equation of linear acoustics, describes an 
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acoustic field resulting from certain specification of the boundary conditions or initial conditions. 
It has little to do with the vortical motions associated with the fluid turbulence that are the source 
of the propagating wave field. Following Ribner (1962) the fluid pressure is decomposed into its 
convective and propagating parts p = p c + p p where p c satisfies 

- Pc, jj = ( pUiUj ),i } (3) 


and therefore p p satisfies 

c^P P ,tt - p p ,j j = -c^Pctt- (4) 

In the adiabatic limit the right hand side can be written as (—/»,«) and is therefore related to 
changes in volume of the fluid element - the dilatation - that generates the propagating pressure. 
In the region of the fluid turbulence, for low turbulent Mach number, the pseudo-pressure is larger 
than the propagating pressure whose source is from fluid motions. Far from the turbulent portion of 
the fluid, the propagating pressure is the major portion of the pressure field as the pseudo-pressure, 
being associated primarily with the convective motions, decays rapidly. Thus, there is an inner 
region of scale t in which the major portion of the pressure is associated with the vortical motions 
and an outer region, or an acoustic mantle of scale A, in which the the propagating pressure is 
the major component of the pressure field. In the inner region of scale l << A, the sound speed 
is effectively infinite: on a time scale of the flow, signals are felt throughout the region of scale l 
effectively simultaneously. Which is to say that, in the near field, p v satisfies the following Poisson 
equation 

- Ppijj = -c^Pett- (5) 

These ideas, well-known in studies of sound generation, were first understood as a singular per- 
turbation by Landau and Lipschitz (1958). A more formal presentation of these ideas is given in 
the following section. To obtain representations for covariances with the dilatation, only the inner 
solution of the singular perturbation problem, were the pseudo-pressure dominates, is used. This is 
consistent with the observations of Sarkar (1992), Blaisdell and Sarkar (1993), in which it has been 
found, numerically, that the incompressible portion of the pressure makes the largest contribution 
to the pressure-dilatation covariance. 

The third idea is that equations should uniformly approach their incompressible form as the Mach 
number goes to zero with bounded derivatives. These facts are used to produce the gauge functions 
in a perturbation expansion in which the small parameter is related to the Mach number of the 
velocity fluctuations, the turbulent Mach number, M t . This does not limit the theory to low mean 



flow Mach numbers. In general, the aerodynamic problem requires an assessment of the effects of 
compressibility on a flow of arbitrary mean Mach number as felt through the_compressible nature of 
the low Mach number turbulent fluctuations. In this way the turbulent Mach number dependence 
of the covariances is obtained by a systematic and consistent balance of terms in the compressible 
Navier Stokes equations. These are ideas have received additional amplification in the very useful 
and thought provoking work of Zank and Matthaeus (1991). 

The use of the inner solution is a useful approximation in mediums that are finite or infinite in 
extent for covariances involving at least one fluctuating quantity which does not propagate - whose 
source is local. Contributions to the covariances from regions outside of the correlation length, the 
outer solution, are negligible. This is because there is no correlation between the local flow field in 
the region £ with quantities outside the region l. Such is the case for < pd > in which the major 
contribution to the fluctuating p is the local p c . 

This is not the case for covariances of the propagating field such as, for example, the variance 
< dd >, whose far field component may be larger than its near field if the size of the domain, 
D, of the flow D/X >> 1 as is implicit in compressible homogeneous simulations. This raises 
some interesting and subtle ideas related to the physics of homogeneous numerical simulations: 
ideas relevant and important to the interpretation and use of homogeneous compressible DNS to 
calibrate models for flows of engineering interest occurring in finite domains. 

The present article makes use of, what is called, the compact flow assumption: the size of the 
turbulent field, D , is small or on the order of the acoustic scale, D/X < 1. Sound traveling through 
a flow' on scales comparable to the wavelength of the emitted sound will begin to be scattered by 
the vorticity. In addition its accumulated effects on the flow' will begin to modify the flow through 
which it is traveling. It is for this reason that acoustic analogies sometimes fail when they are used 
to predict the far field of an acoustic source after the signal has traversed the fluctuating medium 
for more than a few wavelengths. The present interest is in predicting covariances and this compact 
flow assumption is much less of a restriction than in acoustics as the length scale of the correlation 
naturally filters out signals coming from portions of the domain that are uncorrelated with the local 
vorticity. The compact flow assumption is still necessary in order to neglect covariances between 
propagating fields, such as < dd >. 

While most flows of engineering interest can be categorized as compact flow problems, homogeneous 
compressible simulations do not fall into that category. Homogeneous DNS correspond, locally, to 
a turbulence immersed in a general random background wave field which will make contributions 
to the variances of a propagating field such as < dd >, even though the coherence between the 
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local turbulent field and the background field that propagates through it is small. For that matter, 
the local flow could be exactly incompressible and it would still experience, a net drain of energy 
from the background dilatational wave field in which it is immersed. The work of Sarkar et al 
( 199b) has made progress on problems of this type. The model problem they appear to have solved 
in their approach is that of a turbulence of scale l irradiated by an infinite external acoustic field 
generated by a turbulence whose statistics are the same as those of the local turbulent region. 
Their findings are related to flows in which the size of the turbulent field is, at least, A and the 
field is also homogeneous on the A scale (at least) because the assumption of homogeneity made in 
compressible DNS. While such simulations shed much basic insight on the effects of compressibility 
(this paper could not have been written without those insights) application of such results to models 
for compact flows of general engineering interest should be done with these possibilities in mind. 

The assumption of homogeneity is made throughout the mathematical development: this is an 
assumption of homogeneity on the scale l which is to say IjL < 1 where L is the scale of the region 
of homogeneity. This assumption coupled with the compact source and compact flow assumptions 
means that the theory is applicable to flows in which £/ A < 1,Z/A < 1,D/A < L The article of 
Sarkar et al. (1991a) and compressible homogeneous simulations in general appears to be relevant 
to the problem in which l/\ < 1, and L/X> 1,D/A > 1. 

This completes an intuitive background of the physics of the problem and is a useful perspective 
from which to view subsequent developments and more subtle side issues. A mathematically more 
formal statement of these ideas is now carried out. 

2.2 The governing equations: a mathematical foundation 

The following equations are used to describe the portion of the flow of interest : 


Pit T UpPip — P^pip 

(6) 

P^iit pUp^iip + Pii — 9 

(7) 

p/Poo = (p/pccT 

(8) 


For clarity of exposition the viscous terms are not carried: they can be shown to be of higher order 
for the compressible portions of the field, see for example Zank and Matthaeus (1991). This reflects 
the fact that the inner solution of the sound generation problem, on the small length scale, £, is 
being sought and at these scales viscous effects which attenuate wave propagation are unimportant. 
Moreover, a spectral Mach number exhibits an approximate k - 1 / 3 dependence and the scales of 
the motion responsible for the fluctuating dilatation will not be the scales of the flow influenced by 
viscosity. 
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The momentum and continuity equations can be combined to give the following equation 

P,tt ~ P,jj= (pUiUj),ij ( 9 ) 

which becomes a wave equation for p or p if the gas law is used to eliminate one in favor of the 
other. The equation is left in this form for subsequent purposes. There are, of course, some 
limitations regarding the application of this set of equations to a general compressible flow. The 
most substantial is the assumed form of the gas law valid for isentropic flows: heated wall-bounded 
flows have strong temperature effects and as such the application of the theory in the near-wall 
region at high Mach number is, strictly speaking, not valid. The development presented here 
suggests a procedure for handling the problem in more complex flows with heat transfer. 

Perturbing about a quiescent state, (poo,Poo ), the nondimensional forms of the pressure and density 
are taken as p = p^l + p'), p = /^(l + p'). After rescaling the independent variables with t/u 


and t. and dropping the primes, the equations become 

P,t+u p p, p = -(l + p)u p ,p (10) 

(1 + /»)«*»< + (l + p)u p Ui, p +e~ 2 p» = 0 (11) 

P~1P = l/ 2 7(7-l )p 2 (12) 

p, tt -r 2 p,jj = [(1 + p)uiUj],ij (13) 


where e 2 = 7 Mf and M t = u/c where u = 2k/ 3 =< UjUj > /3 and = 'yp <x /p rx> . Note that 
the choice of time scales is determined by the energy containing scales of the motion: it is a coarse 
grained time scale. The fine grained time scale of the problem includes some very interesting physics, 
but relevant only in a cumulative, to the construction of a statistical model. A meaningful balance, 
giving bounded first derivatives of the velocity, is established if p ~ e 2 . It then follows that p ~ e 2 
also. The conventional definition, in concordance with the acoustics literature from which some of 
our ideas are drawn, of the Mach number is used. It is the small parameter that emerges naturally 
in the relevant nondimensionalization of the compressible equations. The use of this symbol as the 
small parameter is only for this section; it will be used subsequently to denote the dissipation. Note 
that the conventional definition of the turbulent Mach number M t = (2Jb/3) 1 / 2 /c means that it is a 
factor 0.816 or 0.577 smaller than the Mach numbers defined using k or q 2 =< ujUj >. Expansions 
of the form 


P = e 2 [pi+( 2 P2 + ■■■■ 

P = e 2 [ Pi + e 2 p2 + •••• 

Ui = Vi + e 2 [ tv, + e 2 w 2 i + .... 
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(14) 

(15) 

(16) 


are chosen. The gauge functions for the velocity are determined by the boundedness condition on 
the second derivative of the velocity, following the methodology of Zank and Matthaeus (1991), 
with different results. The only meaningful balance using a perturbation series in unit powers 
produces the same results. Note that this is not the linear acoustic scaling in which pressure and 
density disturbance scale as p' ~ PooU 2 ~ Pocw'cqq , p f ~ PoqM 2 and thus the fluid velocity of a 
fluid particle associated with the passage of a wave is w f ~ Mu. 

Inserting the expansions into the equations produces, to the lowest two orders, the incompressible 
equations 

Vi it T VpV{,p T pi — 0 (10 

Pujj = - (18) 

Pi = 7 Pi (19) 

where i/j, t - = 0. Note that if the pressure fluctuations are scaled with velocity fluctuations then the 
last equations can be written as pi — M 2 p\. The correction for the compressibility of the flow, 
which does not involve a wave equation on the inner scales, is 


Pi it T Vppi ,p 

= -Wk,k 

(20) 

H T VplViip -\-WpVi ,p T P2ii 

— Pi ( it "t" Vp i ip ) 

(21) 

■ <-> 

ft. 

1 

= ( WiVj + WjVi + Pi ViVj ),ij -pi ,tt 

(22) 

P2 “ 1P2 

= 1/2 7(7 -lR 

(23) 


This is a statement of the fact that, over a region of size £, the pressure signal is felt, effectively, in- 
stantaneously. Reflect on the fact that e is a small parameter and that analysis will not adequately 
represent the effects of compressibility when shocklets are important. The analysis represents the 
effects of compressibility as a linear correction to the nonlinear zeroeth-order of incompressible tur- 
bulence problem. This completes the derivation of the evolution equations for the inner expansion. 
The full problem is the sound generation of acoustics and it requires matching the inner solution to 
an outer solution. For the single-point turbulence closure problem for the dilatational covariances 
the outer solution is not required. The above equations will be used to develop representations for 
the unknown terms. Some of the above relations may also be used to specify initial conditions for 
DNS that seek to investigate compressibility effects whose source is the turbulent velocity field and 
not initial conditions that reflect some other generation mechanism (for example, passage through 
a shock). 

The terms that are sought in this study are various moments of the fluctuating dilatation d = 
Wj,j. The zeroeth-order equations show that the density fluctuations are given by the pressure 


11 



fluctuations, 7 = pi. The evolution equation for the density fluctuations now becomes a diagnostic 
relation for fluctuating dilatation, 


~7 d = P,t+v p p, p . (24) 

The subscript has been dropped. It is seen that one does not need to obtain a solution to the 
evolution equation for the compressible velocity field, Wi. in order to obtain its dilatation. A very 
nice result indeed it forms the kernel of the present pseudo-sound theory. The dilatation is diag- 
nostically related to the local fluctuations of the pressure and velocity; it is the rate of change of the 
incompressible pressure field Pi,jj= («i»j),ij> following a fluid particle. The pressure fluctuations 
which originate as a constraint to keep the eddies incompressible drives the near field dilatation. 
Note that scaling the fluctuating pressure with the mean energy of the velocity fluctuations indicates 
d- M? [p,t+v p p.p ]. 

Constitutive relations for the pressure-dilatation and the dilatational squared covariances can be 
written by taking the appropriate moment of the fluctuating dilatation equation to produce, drop- 
ping the subscript. 


-27 <pd > = < pp >,, (25) 

7 2 < dd > = < pp > + 2 < pv g p, q > + < v p p, p v q p, q > (26) 

for a homogeneous turbulence. The overdot is used to represent the time derivative when it appears 
within the brackets. It should be noted that the near field compressibility effects, as manifested in 
< pd > and < dd >, have been directly linked to the incompressible velocity fields. This fact will 
be exploited to obtain expressions for the dilatational covariances in a turbulent flow. 

3.1 Analysis for the pressure dilatation covariance in isotropic turbulence 

The simplest form of the problem is now solved: expressions for the dilatational covariances for 
an isotropic turbulence without any mean deformation are obtained. For those familiar with the 

pressure strain covariance modeling in incompressible turbulence, this is analogous to the slow 

pressure component of the representation. The constitutive relationship for < pd > is the starting 
point. Using the assumptions of isotropy, homogeneity and quasi-normality, an expression with no 
undefined constants can be obtained. The methods of statistical fluid mechanics similar to those 
of Batchelor (1951, 1953), Proudman (1952), and especially works of von Karman and Howarth 
(1938) are used. 

f A similar expression in different contexts with different assumptions has been obtained independently by both 
S. Girimaji (1995) and S.Crow (1970). 
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Batchelor (1951) has obtained a representation for the pressure variance, < pp >, in an isotropic 
incompressible turbulence. Here, a simpler Greens function method, following Kraichnan (1956), 
is used. The pressure of interest satisfies the Poisson equation: p{xA), 3 j— -(Ti'tj)wj from which it 
follows that the two-point pressure variance obeys 

p{ X . t )p{ X , £) 2 > ijjp'p' — ^ ^iVjVpVq y > ^ijp , q r ( 2 /) 

which can be written as a differential equation in r t = x\ - x x . Following the usual methods for 
translationally invariant random processes, 


< p(x,t)p(x',t) >,jjpp - < ViVjVpVg >,i Jpq (28) 

Using the assumption of translational invariance (homogeneity) where now < p{x,t)p(x', t) >=< 
pp’ > (r). The Greens function for the equation is -^ | r' - r | and the solution is expressed as 

< PP' >= — / < ViVjv' p v' g >, ijpq | r' - r | dV (29) 

The quasi-normal assumption is used to relate the fourth-order moment to the second-order mo- 
ments. The adequacy of the quasi-normal assumption have been investigated over several years. 
Batchelor (1951) has presented evidence of its adequacy when invoked with respect to the large 
scales of the flow. A spectral version of this assumption is used in the EDQNM theory which since 
its inception, as presented in Orzag (1970), has produced very useful results. The adequacy of 
the assumption for the large scales of the flow has been documented in the several experimental 
works. McComb (1990) gives a summary of these results. In compressible flows the adequacy of 
the quasi-normal assumption for the large scales has been investigated in Sarkar et al. (1991a) and 
Blaisdell et al. (1993). The assumption produces 


< V,Vjv' p v' q > = < ViVj >< VpVq > + < Viv'p >< VjV q > + < V x v' q >< V p Vj > . (30) 

The definition for the correlation, < u;u' >= ^R l3 {r), where k = \ < VjVj >, is used to obtain 

2k 

< Vi VjVpV q itjpq 2 <C V{Vp > ,jq < Vj V q >,ip 2 ( ) Rip.jqRjq.ip- (31) 

The pressure variance becomes 


< PP' >= -2 ( ¥) 2 g 7 / R,p,jqRjq,i P I r' - r | d 3 r’. 


(32) 


Continuity, the fact that Rij,j(r) = 0, has been used. For isotropic turbulence the integral can 
be written in terms of the longitudinal correlation function. ^,From the inceptional paper of von 
Karman and Howarth (1938), the longitudinal correlation, < tq(0)i’i(r) > = < iqi'i > fi r ) = 
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7 f/(r) = ^£ 11 , allows the general two-point correlation to be written as R tJ = - r -£-f + (f + 
±rf')6 t j. The integrand can be written in terms of the scalar function /(£)^ Following Batchelor’s 
(1951) development the fourth-order two-point correlation can be expressed as 

= 2 [2 f" 2 + 2/7'" + ~{'f" + 1/7 

= 2 (33) 

Inserting into the integrand and applying integration by parts successively produces, returning to 
dimensional variables, the Batchelor (1951) result: 

9k r°° e 

< PP > = 2 (j) 2 pl jf U'Hm = (34) 

where 

J r oc 

1 tf 2 (Odt. (35) 

0 

Here, and henceforth, £ = r/l is the nondimensional spatial coordinate such that f f{£)d£ = 1. 
Inserting the result into the constitutive relation for the pressure dilatation, 

-27 < pd> = ^ < pp > (36) 

produces, in dimensional quantities, the following expression for the slow pressure dilatation, 

<pd> s = - \l[ §- t i<P> M?k) (37) 

after accounting for the normalization employed. Here, the undisturbed density and pressures have 
been replaced by by the local mean density and pressures. The turbulent Mach number is defined 
as M f = 7 f/c 2 where c 2 = 7 <p> / P is the local mean speed of sound. 

3.2 Analysis for the variance of the dilatation in isotropic turbulence 

The quasi-normal form of the constitutive relation for the variance of the dilatation is 

7 2 <dd>= < pp > + < v p p, p v q p, q > . (38) 

Starting, once again, from the nondimensional Poisson equation for zeroeth-order pressure field, 
p(x = (i’,t7j),y, an equation similar to the two-point variance of the pressure derived above 

can be obtained for the variance of the time derivative of the pressure by differentiation: 

< p(x,t)p(x ,t) >,jjpp = < (ViVj), t (v'pV' q ),t>,ijpq . (39) 

The assumption of homogeneity has been used and the equation is written in terms of the usual 
spatial difference coordinate, r t . Expanding the products of the time derivatives produces, 

< (ViWj), t{VpV' q ), t >, ijvq = 4 < v t vjv' q v' p >, i]vq (40) 
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and the differential equation for the variance becomes 

< p(x,t)p(x',t) >,jjpp = 4 < VjV p ' > vq < Vjv\ >„ P T (41) 

The fact that < Vjv' q >— 0 for homogeneous isotropic turbulence, as can be seen by using the 
Navier Stokes equations to rewrite vj , has been used. The tensor < i ?,•{?'• > can be written in 
terms of the correlation function, < i \v f - > = < vv > R{j which can be rewritten in terms of the 
longitudinal correlation, /i, where as usual, R t j - -^ 7-/1 + (/1 + \rf[)6ij to produce < v z v[ > — < 
vv > [3/i + rf[\ ~< vv > r 2 (r 3 fi ) f . The bi-harmonic equation for the variance of the time 
derivative of the pressure becomes 

< PP > ijjpp = 8 f < t't > jr[$s(rf'f[ )']'• (42) 

Using the Greens function method and integrating by parts produces, in dimensional form, 

< pp >= 4 < vv > tf'fl (43) 

An expression for the two-point variance of the acceleration, < vv > f[ is required. Little is known 
of the longitudinal correlation for the acceleration. The Navier Stokes equations can be used 
to obtain an equation relating the acceleration correlation, / 1 , to /, the well-known longitudinal 
correlation of the two-point velocity correlation K The dynamical equations of the inviscid portions 
of the motion, in the absence of a mean velocity field, can be used to produce the following equation 
for the two-point covariance of the acceleration: 

< Viif'i >= -P~ 2 < pp' >.JJ - < iWjv'iv'b >, jk (44) 

w T here the usual nomenclature of the homogeneous turbulence is in effect and the independent 
variable is the two-point separation r z . 


The quasi-normal expression for the last term on the right in the evolution equation for the tw r o- 
point acceleration produces 


< ViVjVi v' k >,jk = 


[< ViVj >< v'v' k > + < ViVi >< Vjv' k > + < v t v' k >< v'v J 

if ,2 P i/" + ;//')]' (45) 


after substituting in terms of the longitudinal correlation. For operations on functions of r, the 
Laplacian can be written < pp f > V7 = r ~ 2 jp( r2 j^ < PP ' >) an d the dynamical equation for the 
tw r o-point acceleration becomes, after one integration, 

< vv > h = --T-4- <pp’> - - If 2 + -//')• (46) 

Poc r “ r d z r 

^This development was indicated to me by Y. Zhou. 
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^From the expression for the two-point covariances for the pressure, Batchelor (1951), the following 
expression can be derived 

Ifr K PP ' >= - 4pl( T? F ? f ' 2dr ' (47) 

Inserting Batchelors expression in the equation for f\ and taking the derivative produces the quan- 
tity required, 

<vv> f[ = - & 2 (ff" + -//" + - \ff')- (48) 

6 r r r L 


Inserting the expression for < vv > f[ into the variance, < pp > = 4 p^^f < vv > / 0 °° £/'(£ )f[(£)d£ 
produces, after some manipulations, 


< VV >= 4 O 2 ( — l 3 — IZ = O 2 _ A „2 2 TS 

^ rP ^ ^ P oo l i |2 * 3 / k 2 * 2 ~ ^2 ^ ^ 2 ? 


(49) 


where 


n = / 0 °° + in " + irr - £/m • 


(50) 


In the above expression the usual and empirically verified scaling e ~ fi 3 /7 has been used. The 
characteristic velocity fluctuation will be taken to be u 2 = |A: in which case c = a(-£) 3 / 2 /£. Note 
that the coefficient of proportionality is twice as large as that when u 2 —< U\U\ > is used. The 
integral length scale is identified with the longitudinal correlation, l — L n . 


Work by Sreenivasan (1984) has indicated the utility of this expression for turbulent flows. A 
more recent and very timely article, Sreenivasan (1995), assesses the accuracy of this expression in 
several canonical (incompressible) simple shear flows. For homogeneous shear the data indicates 
a ~ 1-2. For the log layer or wake flows a ~ 4. For flows with smaller microscale Reynolds 
numbers Sreenivasan (1984, 1994) shows that a ~ R ^ 1 . There is also a weak dependence on 
nondimensional shear rate. 


The fourth-order moments in the constitutive expression 7 2 < dd >=< pp > + < v p p, p v q p, q > are 
now treated. Beginning with the two-point statistic and writing it as a function of the separation 
distance, r;, 

< v p p, p v' q p', q > >= -[< v k v q >< pp' > + < v p p >< v' q p' > + < v p p >< v' q p >}, pq 

= - < V k v' q >< pp f >, p , (51) 

where continuity, < v p v' q >, p = 0, and the fact that any isotropic vector is zero have been used. 
Further manipulations and setting r = 0 produces 

2k 

^ ^kVq X pp > y P q — ~ pp ]> ,pp 2kb pq < pp > ipq (52) 
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where b pq is the anisotropy tensor, b z j —< Vj > /2k - 1 <5,y . It is zero in an isotropic turbulence. 
A theory including the contribution of the anisotropy of the turbulence w r ouid be possible assuming 
otherwise. Expressing the two-point covariance in terms of its longitudinal correlation function and 
performing the appropriate differentiations of < pp' >=< pp > P(r) produces for the fourth-order 
moment: 

2k 

-y < PP >, PP = -y < PP > 3-Po (53) 

The second derivative of Batchelors solution for the two-point pressure variance can be used to 
show that Pq = where 

roc i 

H = / 7/^ ( 54 ) 

Jo 5 


The fourth-order moment can be written as 


< v p p, p v q p, q > = 


*£<&>$% = MpUf ) 2 Pm 


54 n 2 

^P 00 


(f ) 2 


e rs rs 

TPh 1 3* 


( 55 ) 


Batchelors result for the pressure variance and the empirically validated scabng c = a(~ ) z i 2 ji have 
been used. The particular form of the expression is chosen in anticipation of later manipulations. 


The above results are substituted into the constitutive equation for the variance of the dilatation 
-7 2 < dd >= < pp > + < v p p rp v q p, q > - to obtain the following simple expression for the slow 
portion of the representation for the variance of the dilatation 


< dd > s - 


£ M ( 4 (£) 2 [/f + 6/f/J]. 


( 56 ) 


The variance of the dilatation scales with the time scale of the large eddies of the spectral cascade. 
The integrals are. typically, order one quantities. Zhou (1995) has determined their value from 
high Reynolds number wind tunnel data. Values are given in the appendix. Following present 
conventions, in which the compressible dissipation is defined as < p > e c = | < fi X dd >, the 
model can be put in a form more pertinent to the evolution of the kinetic energy of the turbulence. 
Using the definition of the Reynolds numbers, Rt = to eliminate the viscosity the model can 

be written in a form appropriate for the kinetic energy equation 

16 M? 




c [I s 2 + 6/j/|]. 


3 q 2 Rt 

A complete summary of the models is given in a subsequent section and in the appendix. 


(57 
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3.3 Analysis for the pressure dilatation subject to mean velocity gradients 

In this section an expression for the pressure dilatational covariances for a. general homogeneous 
mean velocity gradient with no mean dilatation is derived. A similar procedure, in which frequent 
recourse is made to the Poisson equation for pressure, is followed. The Poisson equation now 
involves the mean velocity gradient. The nomenclature used in the pressure-velocity covariance 
modeling in incompressible second-order closures will be followed: the covariance will be called the 
rapid component of the pressure-dilatation covariances. 

The constitutive relation for the pressure-dilatation is 

, D 

-27 <pd> = — < pp > (58) 

The velocity field is partitioned according to the Reynolds decomposition + u t ; the upper 
case denoting a steady mean velocity field with constant gradients the lower case will continue 
to indicate the fluctuating field. The mean strain and rotation tensors are 5 , 7 = \ [K,j+V), t ], 
I Vij — | [I’t ,j — V},;]; W 2 and S 2 denote the traces of the squares of these matrices. The nondimen- 
sional form of the Poisson equation for pressure is p(x,t),jj= The fluctuating pressure 

is know given by the following Poisson equation: 


{ v t \ j T V , Vj + Vi Vj ) . 

Multiplying this equation by a similar Poisson equation for p(x',t) and averaging produces 

< p(x,t)p(x\t) = 4 Vuj V P v < v jn r', p . > + < v x v 3 v' p v' q >, upV . 

Expressing the differential equation in terms of the spatial separation, r t , produces a biharmonic 
equation for the t.wo-point pressure variance 


< pp > ijjqq — 4V 2 ',j ,g < Vj Vq > 


(59) 


The fourth order moment which represents the slow pressure contribution obtained in a previous 
section have been dropped. The Greens function method produces the following solution 


< pp > (r) 


4V V 

V 1 piq 


i / 

8 tt J 


< VjV q 


> 


ftp | 


r f \d' 


> 3 / 

r 


4V V I 

^ r HJ v p iq 1 jqip 


(r). 


(60) 


The dependence of the solution on the mean flow gradients has been expressed. The pressure 
variance is known once a representation for the integral I jqip is found. For a class of turbulent flows 
a tensor polynomial in the anisotropy tensor is a suitable approximation for Ij qip . Ristorcelli et al. 
( 1994, 1995a) includes a discussion of issues related to this assumption. Here, only the zeroeth-order 
term in such a polynomial will be retained for the purpose understanding the physics and obtaining 
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scalings for the compressibility effects. Higher order terms which scale with the anisotropy of the 
turbulence. b{ 3 —< v 3 > /2k — ^<5- are neglected in these zeroeth-order expressions. Note that 
b \ 2 ~ 0.16 is not atypical. At this point, and considering the purpose of the article, the additional 
algebra necessary to obtain the anisotropic contributions to the integral is not warranted. A 
fourth-order isotropic tensor possessing the proper symmetry and satisfying the continuity relation 
Ijiip = 0 at r = 0 is 

Ijqip ~ ^1 \^jq^ip "h ^jp^zq)} 


where 


jsj _ A/ AA f <vv'> r'd 3 r' - — —PF 

1 ~ 15 “15 8 ir 1 J 9 >,tp 15 3 1- 


(61) 


Expressing the integrand in terms of the longitudinal correlation in the normalized coordinate. 
< VjVj >,u = f f _2 [f 3 /'" + 7 £ 2 /" + 8 if]. The facts < VjVj > = f [r/' + 3/] = r~ 2 £{r 3 f) and. 
that in spherical coordinates, the Laplacian is V 2 = r ~ 2 ^ r2 jp have been used. It is also possible 
to integrate by parts allowing the integrand to be expressed in lower order derivatives for more 
accurate computation from experimental data. Thus 


a\ = if e / 0 °° e\er + ur + sm 

TS$*f“Udt=±ft 2 r 1 (62) 


and the solution for the rapid pressure variance in an arbitrary three-dimensional mean velocity 
gradient can be expressed as 

< PP > r = yi 2 ii 2 [35 2 + 5 IT 2 ] /[. (63) 

Note that the integral has dimensions of a characteristic correlation area: the rapid pressure contri- 
bution to the pressure variance will vary according to the spatial scale of the turbulence unlike the 
slow pressure contribution given. This dependence on the spatial scale was first noted by Kraichnan 
(1956) who solved the problem of the pressure fluctuations in an isotropic turbulence in unidirec- 
tional shear. The results here extend Kraichnans (1956) results for a planar unidirectional shear 
to an arbitrary mean deformation. A very interesting, and more modern, paper highlighting the 
physical and wavespace aspects of the results is George et al. (1984). The results are now substi- 
tuted into the constitutive equations given above and the rapid portion of the pressure-dilatation 
covariance in dimensional variables is then given by 

f|V +5H ' 2 i] 
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Note the appearance of the quantity S£/c 0 c ; the dependence of the compressibility effects on a 
deformation rate Mach number indicated by Lele (1994) and Sarkar (1994) and substantiated 
phenomenologically in Sarkar (1994). Note that the theory also predicts a dependence on a Mach 
number based on the mean rotation, W£/c 0 c . The expression is recast in terms of the turbulent 
Mach number and the representation for the rapid component of the pressure-dilatation covariance 
becomes 

< Pd > r = ~ 7 l Wt [< P > [352 + 5W2] ] (64) 

It is seen that two effects contribute to the pressure-dilatation covariance: one due to the exchange 
between potential and kinetic modes of energy (since M 2 ~ k ) and the other due to changes in the 
scale area of the correlation. Increases in the kinetic energy results in a transfer of the mechanical 
energy to the fluctuating pressure field. Similarly, increases in length scale, implying a decreased 
rate of cascade to the smaller scales, also transfers energy to the fluctuating pressure. Unlike the 
slow pressure-dilatation, however, the rapid-pressure-dilatation does not always have the opposite 
sign of the growth of kinetic energy but now depends on the rate of increase of the area of the 
correlation, t 2 . 

In order to close the representation it is necessary to have an expression for l . This introduces 
an element of empiricism; up to this point no phenomenological assumptions, other than the very 
reasonable quasi-normal assumption for the large scales, had been made. The usual heuristic 
approximation e = a(2fc/3 ) 3//2 /f where a ~ 1 — 4 produces 

< Pd > r = ~ ~ (|) 3 /[ o 2 ~ [< p > kM 2 [3S 2 + 5W 2 ] } (65) 

Here the quantities with a carat are nondimensional deformation and rotation rates eg., S 2 = 

(Sk/e) 2 . 

3.4 Analysis for the variance of the dilatation subject to mean velocity gradients 

In the constitutive relationship for the variance of the dilatation the time derivative is replaced 
by the mean advective derivative, ^ = ( ), t + V k ( ), k which comes from the mean portion of 
the advective terms in the expressions for the dilatational covariances. Carrying the substantial 
derivative as part of the time derivative term involves no approximation and follows quite naturally 
from the Reynolds decomposition. However it is necessary to carry out the development in a 
way that preserves Galilean invariance. The quasi-normal form of the constitutive relation for the 
pressure dilatation is 

7 2 < dd >= <pp> + < v p p, p v q p, q > 
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The small circle will be used to indicate the mean convective derivative; for example in the equation 
above <pp>- < 75? >♦ 


An expression for <pp> appearing in the the constitutive relations for the variance of the dilatation 
is obtained first. For the convenience of the presentation the two contributions to the variance of 
the dilatation will be denoted 7 2 < dd and 7 2 < dd > 2 = < v v p, p v q p, q >. Applying 

the Reynolds decomposition to the nondimensional form of the Poisson equation for pressure. 
p{xj),jj= — {viVj),ij, and taking the appropriate derivatives and dropping the terms quadratic in 
the fluctuating velocities (which were treated in an earlier section) produces 

p ( *£ • 0 'fj j ~ ( lj T 1" i — 2\i,j Vj,i • (66) 


Multiplying this by a similar Poisson equation for p (x\t) and averaging produces the biharmonic 
equation for the two-point pressure variance 


<P (.T, t) p (x ,/) ^yj'j'qq — 4 Vj,j Vpiq' Vgip 1 ^ — 4lj, j \ p,g < Vj V ^ > , 


ip • 


(67) 


The last equality has been written in terms of the separation variable, r x . The Greens function 
solution procedure produces the following representation for the two-point variance 

<pp > = 4 Vi,j V p , q J <VjV q >, ip | r - r'\d 3 r' = 41;, j V p , q I jqtp {r). (68) 

It is the variance that is required. Following the method discussed in the the previous section, the 
fourth order tensor, neglecting higher order corrections for anisotropy, is represented at r = 0, as 
an isotropic tensor 


— Ay\6jq$%p . (fiji 6qp 4" ^jp^iq)\ 


l JQ l P --zi-jy-iy ^ 


e 2 r 
15 Jo 


° ° a! 1/ 

^ V j ^ j ^ ’PP £ = 


1 2 2k 


(69) 


15 Jo 15 3 

As has been noted little appears to be know about the two-point statistics of the acceleration. If the 
acceleration correlation were known it would be a simple matter to show that the integrand is given 
by < VjVj >,u =<vv> [£ 3 fi' f + /{' + 8£f[}t- 2 in the normalized coordinate. Unfortunately this is 

not the case and an expression for }\ in terms of / is sought . The Navier Stokes equations without 
the viscous terms, which describe the energy containing range of the flow, will be used to obtain an 
expression for the integral, <VjVj>, pp { 3 d£. Taking the equation for - ij = VkV{,k +VkVi,k+Pi, 

o f t 

and multiplying it by a similar equation for v -, averaging and taking the trace produces, in the r t 
coordinate, 

<VjVj>= - Vi,k Vi,g < VkV'g > - < pp' >,jj - [V;,fc < Vkv' q v\ >,, +Vi, q < V k V;V q >,k ] 

- [< P v 'k v 'i + < P'ViVk >nk ] - [V;,* < V k v' q v' t > , q + V;, g < V k V,v' q >, k ] 

- < Vj VkV q Vj >,kq 
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The two-point triple covariances are zero for homogeneous isotropic turbulence and the fourth-order 
correlation was treated in a previous section. The equation yields, after taking the Laplacian, the 
quantity sought: 


o o 1 

< v j v j>,pp — 


^itk. i hij ^ Vq ^ ipp ^ pp > ijjpp * 


(70) 


In the previous section it had been shown that the two-point pressure variance satisfied the bihar 
monic equation: < pp f >,jj qq — 4 V^j V v , q < VjV f q > , tp . Thus 

^ i ik ^ Vq ipp 41';', j I'p.g K Vj V q > ,jp , 


<VjVj> ipp 


71 ) 


which upon multiplication by £ 3 and integration produces the desired result for / 0 °° < vj v- > , pp £ 3 d( > 


in the definition of /£: 


2k 

~~z~ * 0 


^iik iiq ^kq 4 V t ' , j V p , g I 


j q ip 


(72) 


and I\ is seen to be related to the two integrals, I kq = / 0 °° < v k v f q >, pp £ 3 d£, and I jqtp = / 0 °° < 
VjV f g >,i p t; 3 d£. The isotropic portion of these tensors are related to an earlier integral, /[, defined 
in the previous section. The tensors have the following representations 

2 1 

I jqip “ 15^1 4" fijpfiiq)] 

l jq — gM V 3V 

Inserting these expressions into the equation for L \ produces an expression for the two-point accel- 
eration correlation integral in terms of the two-point velocity correlation integrals: 

Ci — ^ [135 2 + 15W 2 ] , (73) 

and the rapid pressure variance becomes 

<°PP >r = TT [35 2 + 5W 2 } [13S 2 + 15 W 2 } I\t 2 y . (74) 

Substituting l = a(2k/3) 3/2 /e, and inserting into 7 2 < dd >i=<pp> which is related to the 
dilatational dissipation by 4 = & < dd >! = <pp> produces 

1 9 A/f 4 

4 = (^(j) 5 ^ e [3^ 2 + 5^ 2 ] [1.35 2 + 15W /2 ] a 2 I[ (75) 

after accounting for the nondimensionalizations employed. 


An expression for the fourth-order moment, < v p p, p v q p, q >, appearing in < dd > 2 is now sought. 
In a previous section it was seen that under the quasi-normal and isotropic approximations that 
< v p p, p v q p, q >= — ^ < pp >, qq . The Greens function method produces 


, -If , d 3 r' 

< PP > vj ~ 41 V? 47 _/ < v i v q > »p i r _ ~i 


|r — r'| 


\ p , q Ijqip( r ) 


(76) 
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where the biharmonic equation for the pressure variance from the previous section has been used. 
Following the usual procedures, at r = 0, with I jq i p ~ [6j g 6i p - \(6ji6 qp ± 6j p Si q )\ produces 

9 9k 

<PP>,n= -y [3S 2 + 5W 2 }F 3 . (77) 

where 


A 


r 

3 — 



2 2 k 


T r 


15 3 


(78) 


Using the facts < VjVj >= ^[r/' + 3/] = r _2 ^:(r 3 /) and, that in spherical coordinates the 
Laplacian is V 2 = r ~ 2 ^ r2 j^^ produces /£ = - / 0 °° £ 2 / / " + 7£/" + 8/' d£. Reflection on the integrand 
will show that it is suitable for application of Gauss’s theorem: the exact result I JJlt =< VjVj > = 
~3/(0) is possible and 1 3 = 3. The integral is nonetheless carried symbolically in the light 
of computations extending this theory to anisotropic turbulence. More pragmatically, the integral 
expression and its exact value, have been used to evaluate the accuracy of the numerical integration 
technique of experimental data. 


The fourth-order moment becomes 

2k 2 2k 

< v p p, p v q p, q >= - — < pp >,„ = — ( — ) 2 [3S 2 + 5 W 2 ]I$ (79) 

6 1 5 6 

and thus 


< dd > 2 = [35 2 + 5W /2 ]/ 3 r . 


(80) 


which allows the second portion of the rapid dilatational dissipation to be expressed as 


c c2 




(81) 


The rapid portion of the dilatation dissipation can be written as the sum e r c = e T cl + e r c2 and thus 
9 M 4 q 1 

< = (|) 5 ^ [3^ 2 + + (^) 2 [13S 2 + 15ir 2 ] a 2 /[]. (82) 

This concludes the analytical development of the representations for the pressure dilatation and 
the dilatational dissipation covariances for compact flows. There will be additional mathematical 
manipulations of the expressions obtained section in order to understand the implications of the 
analysis. 
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3.5 Summary of the dilatational covariance representations 

The section is ended with a summary of the results of the analysis. The dilatational dissipation is 
comprised of a slow' and a rapid part: e c = e r c where 
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(83) 

(84) 


It is useful to reflect on the results of Blaisdell et al. (1993), Figure 12, in which it appeared that 
the dilatational dissipation could not be parameterizable solely in terms of the turbulent Mach 
number. The present analysis suggests its dependence, in simple shear flows, on two additional 
parameters, R t and Sk/e s . The pressure-dilatation covariance is a sum of similar terms, < pd >=< 
pd > s + < pd > 5 , 


< pd > s = ~ \ll T [< p > M?k] (85) 

< Pd > r = - ^ (|) 3 /[ [< P > kM 2 [3S 2 + 5 W 2 } ] (86) 

The constants, denoted by the / n in the these expressions are given by integrals of the longitudinal 
correlation: 


r oo 

n = / u ,2 dt 

Jo 

n = - jH uvr + i/r + pr - ^/m 

« - /o «• r d( 
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I\ = 2 / 
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n = - er+ur+sf di 
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Except for two very reasonable phenomenological assumptions the results presented above are 
a mathematical consequence of the assumptions that led to the diagnostic relationship: —7 d = 
Pit+VpP^p. The assumptions used in developments subsequent to the diagnostic relationship are 
the quasi-normal behavior of the large scales and relationship relating length scale to dissipation. 
The analysis, apart from the quasi-normal behavior, verified in Blaisdell et al, (1993), has produced 
an exact but unclosed, in the context of single-point moment methods, result. The quasi-normal 
approximation relates the fourth-order moments of the velocity distribution to the second-order 
moments as if the large scales of the turbulence were Gaussian. To achieve closure an expression 
for the length scale is required; the very well established phenomenological relationship between 
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turbulence length scale and dissipation: t — a(2k/3) 3 / 2 /e s is used. With these qualifications in 
mind the expressions derived are, in the limit of a homogeneous isotropic turbulence, mathemat- 
ically precise. The expressions here may be viewed as the leading order term in a more general 
expression in which successive terms scale with the anisotropy and inhomogeneity of the flow. It 
is expected that such an analysis for the dilatational covariances will, at the very least, predict 
the fundamental nondimensional parameters and scalings in the characterization of the effects of 
compressibility. 


4.1 The physics embodied in the dilatational dissipation representations 

Various aspects of the compressible dissipation representations are now discussed. The dilatational 
dissipation is comprised of a slow and a rapid part: e c = e r c + e s c where 


16 Mf 
3 ^ "fiT 


[/! + 6 im 


: \f ^ fj [35 2 + 5 W 2 ][\ll + (-V 2 [13S 2 + 151V 2 ] o 2 /[] 
3 Kt 5 la 


(87) 

( 88 ) 


Immediately apparent, in contradistinction to other models for these terms, is the fact that the 
analysis' predicts a dependence of the compressible dissipation on mean flow gradients and Reynolds 
number. The dependence on the Reynolds number suggests that assessing the importance of 
the dilatational dissipation on the basis of low Reynolds number numerical simulations may be 
misleading when applied to higher Reynolds number flows. Computation done with these models, 
for example Ristorcelli et al (1995) indicates that the major reduction in spread rate in the mixing 
layer, for example, is due to the pressure dilatation. The Mach number dependence is also stronger 
than the M 2 dependence in Sarkars model for the dilatational dissipation, and less steep than the 
exponential dependence of Zemans model. 


The importance of the dilatational terms is difficult to assess a priori ; their scaling with M t and 
R t suggests that they are negligible. This is probably the case for the slow term. This is however 
not the case for the rapid term. Terms like S appearing in e r c are typically in the range 0 - 10; 
depending on the mean velocity gradients the dilatational dissipation may or may not contribute 
to a flow. The equilibrium (incompressible) homogeneous shear, for example, has Sk/e ~ 6; the 
equilibrium log-layer Sk/e ~ 3.3. For these planar flows S and W are the same. Using 5 ~ 10 as 
an upper bound the quantity [35 2 + 51V 2 ] ~ 10 3 . When squaring it again, as occurs in the /[ term, 
the fact of its M 4 dependence is easily compensated. 

The present mathematical development also shows the importance, in mean fields solely charac- 
terized by the mean strain and rotation, of a Mach number based on the mean velocity gradient. 
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Sarkar (1994) has defined a gradient Mach number as M g = Sl/c where, for simple planar flows, 
Ui,j — U \,2 = S and the length scale he used is transverse two-point correlation of the longitudinal 
velocity. In this article i will be taken as the traditionally defined integral scale. Using the scal- 
ing l — a(2fc/3) 3 / 2 /e s the dependence on the mean deformation can be related to the turbulent 
Mach number and the ratio of strain to correlation times: Stjc — or| = a| SM t = Ms- 

For a general three-dimensional flow the theory makes a distinction between a Mach number 
based on the mean strain and the mean rotation. A second gradient Mach number is defined 
Wife = a §T T'Mt = ot\WMt — M\v. For simple planar flows the mean rotation and strain 
are the same, M$ = Mw and the mean gradient parameterization, St/c is complete. For flows 
characterized by mean pressure gradients and bulk dilatation additional work is required. 


The rapid portion of the dilatational dissipation can be rewritten in terms of these two mean 
gradient Mach numbers as 


M? 

Rt 


IT" €s + 5Mw][-Iz 
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+ (— ) 2 [135’ 2 + 15W" 2 ] a 2 /[] 
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(89) 


The structure of the model is seen to be similar to the Sarkar model, e c a M 2 e 5 , but now with a 
coefficient that is not a constant but depends on the Reynolds number and the mean strain and 
rotation Mach numbers. 


Sarkars (1994) heuristic reasoning can be used with success to indicate the behavior of the dilata- 
tional terms in different mean flows. Though Sarkars (1994) subject is the changes in the anisotropy 
of the turbulence due to compressibility, as indicated by the work of Abid (1993), his arguments 
are equally applicable to both the dilatational dissipation and the pressure-dilatation covariances. 
Sarkars (1994) arguments indicate that the effects of compressibility are much larger in the mixing 
layer than in the equilibrium boundary layer: the mixing layer is stabilized with respect to the 
boundary a layer by compressibility. The difference between the compressible mixing layer and the 
boundary layer flow can be parameterized in terms of a gradient Mach number, M g . In Sarkars 
(1994) examples M g (proportional to Ms) for the mixing layer can be an order of magnitude larger 
than that for the boundary layer. The same reasoning using the a mean gradient Mach number 
applied to the dilatational dissipation indicates that compressibility dissipation effects are substan- 
tially more important for the mixing layer than for the wall boundary layer. Using Sarkars (1994) 
values and definition of the gradient Mach numbers, M g ~ 6 in a mixing layer while in the boundary 
layer M g ~ 1 and the effects of the compressible dilatation are an order 6 2 more important in the 
mixing layer than in the boundary layer. A similar variation is seen using Ms to characterize the 
effects of compressibility. 
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Gradient Mach number type quantities have been identified in the works of Durbin and Zeman 
(1992), Cambon et al. (1993), and Lele (1994). The quantity has been ,given various physical 
interpretations. Durbin and Zeman (1993) have, in the context of a RDT theory for compressed 
turbulence, interpreted it as the change in Mach number over an eddy of scale l . It can also be 
understood as the ratio of the acoustic propagation time across an eddy to the mean deformation 
time scale, Lele (1994). Reviewing the analysis indicates a more thermodynamic, rather than kine- 
matic, origin and interpretation of the gradient Mach number. The St in the gradient Mach number 
appears because of the scaling of the rapid pressure integral with the area, as was first noticed by 
Kraichnan (1956). The c in Sl/c appears because of the linearization of the gas law relating the 
fluctuating p and p about, the local mean pressure and density state. It sets the magnitude of the 
proportionality constant in the dimensional form of the diagnostic equation relating d and p. With 
these ideas in mind one is led to interpret the gradient Mach number as an indication of the relative 
magnitude of the pressure fluctuations (due to shear ) to the dilatational fluctuations as set by the 
local mean density and pressure (which need not be adiabaticallv related). 

Mention should also be made of the Wilcox (1992) analysis of the sensitivity of the flat plate friction 
coefficient, c/, to models for the dilatational dissipation. His arguments show that models, such as 
the Sarkar et al model (1991b) or Zeman (1990), predict effects of compressibility when in fact 
there are none. Sarkars (1991b) model, for example, undesirably reduces the skin friction in the 
compressible flat plate flow because of the modification of the effective von Karman constant. (It 
should be made clear that Sarkars model was intended for use in free shear flows such as the mixing 
layer: it was after all calibrated using the homogeneous DNS of a compressible shear). Wilcox’s 
analysis has been repeated for the current model. The modifications to the von Karman constant 
of the present model are smaller than that of the incompressible form of the modeled equations. 
This is because of the near wall Mf dependence. 

The thoughtful reader will have noticed that the analysis has produced a representation for the 
dissipation that depends on the Reynolds number. The magnitude of the dilatational dissipation 
depends on the viscosity: for a fixed Af*, as R t — ► oo, the dilatational dissipation vanishes. This is a 
rigorous consequence of the diagnostic relationship, - 70 ? — p, t +v p p , p , derived from the perturbation 
method and subsequently employed to obtain the results. The initial assumptions lead to an 
expression for the dilatation that is related to the pressure fluctuations, an essentially inviscid 
phenomena. The dependence on the viscosity arises when one computes the compressible dilatation 
from its definition using the variance of the dilatation: e c — \v < dd >. 
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It appears that, in the small turbulent Mach number limit, that the usual interpretation of dissipa- 
tion type quantities as spectral fluxes is not appropriate for the dilatations! dissipation. This is a 
finding which has caused some consternation and much consideration; it is, however, a mathemat- 
ical consequence of the initial assumptions. A portion of the ideas arising from diverse discussions 
with colleagues will serve to make this result plausible; as will, perhaps, the fact that the results are 
consistent with EDQNM results. The utility of the results of the analysis, in the context of their 
application to computing engineering flows, is another and very different subject to be treated in a 
subsequent work. 

That the compressible dissipation might not be interpreted as a spectral flux (as is the case of 
the solenoidal dissipation) is suggested by results given in the EDQNM of Bataille (1994), In 
Batailles (1994) simulation the energy spectrum is divided into its incompressible components and 
compressible components. The solenoidal spectrum, E ss , is found to scale, as is usual, k“ 5 / 3 ; the 
compressible spectrum, for small Mach number, is much steeper and scales as E cc ~ k -11 / 3 . 

Multiplying by k 2 the solenoidal and dilatational dissipations are found to scale as k 1//3 and k~ 5 / 3 . 
The negative power law scaling of the dilatational dissipation indicates that, unlike the solenoidal 
dissipation, the dilatational spectrum peaks in the lower wavenumber regions of the spectrum. 

Such a point of view can also be understood by more heuristic arguments involving a spectral Mach 
number, M 2 {n) ~ E(k)k/c 2 . Using the incompressible spectrum, as the compressible spectrum 
falls off faster, produces M t (k) ^ k ~ ly/3 suggesting that the dilatational dissipation is a result of 
a combination or competition of effects that are important at different scales of the motion: the 
energy in fluctuating dilatation at the large scales and the sharp gradients necessary for viscous 
dissipation at the small scales. This is intuitively consistent with the fact that, for fixed M t , 
increasing R t by decreasing the viscosity adds more small scales to the field that are also more 
divergence free. Thus, with the length at which the gradients are strong enough to undergo viscous 
dissipation becoming smaller and simultaneously more divergence free, there results a net reduction 
of the dilatational dissipation. 

These ideas must be tempered with the fact that these results come from a low turbulent Mach 
number perturbation method - a linear perturbation about the nonlinear incompressible problem. 
This is not to say that physics related linearly to the velocity field exhibit a steeper spectral slope 
than the energy spectrum. The linearly related spectrum of a passive scalar temperature, which 
scales with ~ k“ 5 / 3 , Tennekes and Lumley (1972), doesn’t. The problem is more complicated than 
this and the analogy is inappropriate. While a low turbulent Mach number analysis is expected to be 
appropriate for most flows of aerodynamic interest, it must be remembered that the nonlinear self- 
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interaction terms of the compressible velocity field are absent in this analysis. As M t — *• 1 nonlinear 
terms and shocklets are expected to become important. This analysis has no relevance to flows in 
which shocklets are a major portion of the dissipation. On this point the ones attention might go 
to the homogeneous shear results done as part of G. BlaisdeLls thesis, published in Blaisdell et al. 
(1993); it was found that shocklets, for moderate turbulent Mach numbers, contributed very little 
to the dilatational dissipation. Note the difference in definition of the turbulent Mach numbers. 

Batailles (1994) EDQNM results show that as the M t increases the slope of E cc decreases. Though 
outside the range of the validity of the simulation, the EDQNM simulation shows that as M t — *• 1, 
the slope of E cc is approaches the slope of E $s . These speculations concern phenomena outside of 
the range of validity of the linear theory and the EDQMN; within the range of validity the results 
of the perturbation theory are consistent with the EDQNM simulations. These EDQNM results 
have been found to be insensitive to the form of the small scale damping terms. 

4.2 The physics of the pressure-dilatation covariance representations 

The phenomenological implications of the pseudo-sound assumptions for the pressure-dilatation 
covariance are now explored. Unlike the compressible dissipation, which represents an irreversible 
transfer of energy, the pressure-dilatation represents a reversible transfer of energy between kinetic 
and internal modes. This reversible rate of transfer is proportional to the departure of the flow 
from equilibrium and with a simple rearrangement of terms the pressure- dilatation is seen to be 
equivalent to an increase in the flows inertia. Subsequent to this discussion the representations 
for the pressure-dilatation covariance are further manipulated to produce a final expression that is 
more easily understood and applied. The pressure-dilatation is seen to be proportional to the net 
imbalance of production, transport and dissipation of k and T. The pressure dilatation is also seen 
to be a function of the how rapidly the eddy turnover rate, k/e s , tracks then mean deformation 
and rotation, S and IT. 

The pressure-dilatation covariance is a sum of the slow and rapid terms < pd >-< pd > s + < pd > s 
already given. The full pressure-dilatation covariance representation can be written as 

<pd> = -I pd §;{<P>kM?] -r pd R[ <p>k MfT] (90) 

where 


T = 

[35 2 + 5 IT 2 ] 

Ipd 

+ I r pd T 
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Let T be called the relative cascade or eddv turnover rate. By expanding the differentials according 
to the product rule the expression for < pd > can be suggestively rearranged in the k equation to 
produces terms representing an “added mass” effect as well as two additional source terms: 

[1 + IpdMf] < p> k = Pk ~ < p> e - k I pd j^[< p> Mf]- < p> k Mf I r pd T. 

The evolution of the kinetic energy experiences additional inertia. The flow appears to act with an 
added mass equal to Mf times the weighted sum of the integrals: / p j and I^ d . The additional inertia 
is due to the reversible transfer and storage of energy in the internal energy (mean temperature) 
field. This role of the pressure dilatation, as a transfer between internal and kinetic modes of 
energy, appears to have first been noticed by Zeman (1991) and explored further in homogeneous 
shear by Sarkar et al. (1991a). Thus the effect of the pressure-dilatation is to reduce the effects of 
production and dissipation unbalance by a factor 1 + I p dMf. 

0 0 

Note also the appearance of the nonequilibrium and history effect, T and M t . The first reflects how 

O 

rapidly the eddy turnover time tracks the mean velocity gradients. While the second, M t reflects 

O 

how rapidly the kinetic energy and mean internal energy adjust to each other. Thus k is influenced 
by the rate at which the pressure dilatation can equilibrate the “potential” difference between the 
k and T fields. The appearance of relaxational effects incompressible flows have been noted in the 
calculations of Abid et al. (1995). 

The mean deformation and mean rotation rate Mach numbers play a role in the added mass term. 
They appear in the product of I v dMf. Thus, the gradient Mach numbers, as manifested in the 
term MfT, also effects the development of the flow by influencing the pressure-dilatation terms. 
The representation distinguishes the mixing layer and the boundary layer not only through the 
substantial derivative terms but also through the mean gradient Mach numbers. These terms are 
expected to make a difference primarily in flows with streamwise variations. 

Additional analysis will now produce a final form for the pressure-dilatation representation. The fact 
that the pressure-dilatation covariance depends on the rate of change of Mf shows that represen- 
tation couples the kinetic energy equation to the internal energy equation. The pressure-dilatation 
covariances can be expressed in terms of k and T through the definition of Mf; this suggests that 

O 

Mf in the representation can be eliminated. The substantial derivative of Mf is easily found from 
its definition and the ideal gas law to be 

O 0 

M?= Mf [ | | ]. ( 91 ) 
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The mean temperature and k equations are, with the assumption of constant c v , 

< ft > c v jj- t T = Pp-<pd> + <p> c 5 + < p > € c +- Tj 

< p > -jj- t k = P k + < pd > - < p > e s - < p > ( c + T k 

O 

will be used to eliminate M t from the equations. Examination of the two equations for k and T 
shows some important effects. Both T and k receive energy from the mean flows kinetic energy 
through the production terms, Pj and P^. The signs of PD and P ^ can be of either sign in 
which case energy from the mean flow can be diminished or increased by interactions with the 
turbulence or with the mean internal energy of the flow. Note that if one insists on a Boussinesq 
eddy viscosity approximation for the Reynolds stresses that P \ allows only a one way transfer of 
energy from the mean to the turbulence. The dissipation terms, < p > e s -\- < p > e c , are always 
positive representing a flow of energy from the turbulence to the internal energy of the fluid. The 
pressure-dilatation, exchanging energy between kinetic and internal modes, on the other hand can 
be of either sign. 

The pressure-dilatation covariance, for a flow with negligible mean dilatation, can be written as 

o 0 

<pd> = - <P> IpdM? k[ 2^ - ^ ]-</>> k M? Ip d T (92) 

O O 

Using the equations for k and T and the definition of the Mach number in terms of temperature. 

— §A/* 2 t (7 — 1) produces a simpler and final and almost algebraic expression for the pressure- 
dilatation covariance: 


< pd >= -XpdM? [P k - < p> e + T k - 7(7 - 1 )(Pt+ < p > ( + T t )] 

- < p> k M? x r p d T- (93) 


The stream wise adjustment of the eddy turnover time scale to the mean velocity gradients remains. 
Here e stands for the combined solenoidal and compressible contributions to the dissipation and 


Xpd 

Xpd 


2Ipd 

T T 

pd 

l+2/ pd M ( 2 + f/; d M t ‘‘7(7-l) 


The \ coefficients are thus functions of the turbulent Mach number, , and the relative turnover 
rate T through I p d . 


One of the failures of early models for compressibility effects in turbulent flows is their lack of 
universality. Current models capture the compressibility effects in the mixing layer, but the same 
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models, when applied to the near equilibrium boundary layer, add undesirable compressibility 
effects. Here the dependence of the results on the mean flow parameters, such as the mean strain 
and rotation Mach numbers and in DjDt - as manifested in the production-dissipation balance and 

o 

in T - show that the analytical representations do distinguish between these two classes of flows. 

The effects of compressibility have been observed to be negligible in the unidirectional near equi- 
librium boundary layer flow. In near equilibrium flows there is an approximate balance between 
production, dissipation and transport and therefore DjDt ~ 0: the model then predicts, consistent 
with observations, < pd 0. Particularly noteworthy are the predictions for the equilibrium log 
layer in which the transport terms negligible: for the log layer Pk ~< P > € and < pd >~ 0 + Q{Mf) 
is a further indicating that these compressibility effects are not important in equilibrium flows. In 
the mixing layer or jet, on the other hand, where production is an important quantity, the pressure- 
dilatation will shunt energy from the turbulence into the mean temperature thus reducing the level 
of the kinetic energy. This has been seen in laboratory experiments an also in calculations done us- 
ing these analytical results. These calculations are the subject of studies addressing issues relevant 
to turbulence modeling and computations. 

For flows in which production is not important the analysis indicates that the net effect of the 
pressure-dilatation is (if the contributions of the dissipation and transport to the energy budget 
have the same sign) is to increase the level of the the kinetic energy of the turbulence. Such 
situations arise in wake flows with and without momentum defects. 

Expressions for the pressure-dilatation in two simple flows, the isotropic decay and the homoge- 
neous shear are worth considering. Consider first an isotropic decaying turbulence. The pressure- 
dilatation covariance is 

< v d >= x v dMt [1 + 7(7 - 1 )}< p> e. (94) 

It is seen that the pressure-dilatation is positive indicating a net transfer of energy from the internal 
modes to the turbulence. After which the energy is, of course, dissipated by viscosity and returned 
to the internal energy of the fluid increasing its temperature. A portion of this energy, proportional 
to Mf and the extent of the departure from equilibrium, can then once again transferred to the 
turbulence. 

o 0 

The T and k equations for the case of an isotropic decaying turbulence are written 


Dt 

- (1 - XpdM? )c + 

\xpdM} 7(7 - 1)« 
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and the dependence on the production-dissipation balance and the turbulent Mach number are more 
readily seen. Note that the factor multiplying the dissipation is always positive, 1 — XpdM? > 0. 
Here e = e c + e s and that the density has been set to unity for convenience of presentation. Note 
that the energy decay rate is explicitly dependent on the Reynolds number through e c . It is, of 
course, still implicitly dependent on the Reynolds number by the dependence of e s on the Reynolds 
number in a decaying turbulence. This issue has been inconclusively explored in Blaisdell et ai 
(1991) in an attempt to assess the effects of compressibility on the decay law. 

The simulations of the isotropic decay by Blaisdell et ai (1991) and Sarkar et ai (1991) have 
indicated a strong dependence on the initial conditions. Neither of these simulations have used 
initial conditions consistent with the pseudo-sound analysis. This would require initial conditions 
in which 7 < pp >— ~ j < 86 >=< pp > where < pp > is the variance of the incompressible 
pressure field. Which is to say that the compressible fluctuations are generated by the turbulence 
as opposed to imposed on the flow as an arbitrary initial condition. Even the choice of socalled 
“incompressible” initial conditions < 08 > = < pp 0 is asymptotically inconsistent with finite 
initial M t . This can be verified by expressing 7 < pp > = < pp > in primitive variables in which 
pressure is nondimensionalized by 

Blaisdell et al. ( 1993) has investigated the possibility a polytropic gas law, n < pp > = < pp > where 
< pp > is the total pressure variance. It is found that n ~ 7 in the homogeneous shear and that, in 
the isotropic decay, n is initial condition dependent. One may well conjecture that a set of initial 
conditions as specified by this pseudo-sound theory might show n ~ 7 subsequent to intialization. 
In which case Blaisdell et al. (1993) conclusion that algebraic models for the dilatational dissipation 
are inadequate may have to be qualified. We do concur with the Blaisdell et al. (1993) conclusion 
that in situations where the compressible component of the flow is arbitrarily specified by the initial 
conditions that such algebraic models will not work. 

Sarkar et al. (1991) have modeled - < pd > +c c = a\ M?e s and the turbulence energy equation 
can be rewritten 


C p > jyi k =< pd > - < p > e s - < p > e c = -(1 + ai Mf)e s . 

Keeping only order Mf terms the present analysis gives for the turbulence energy equation, in 
apparent contradiction 

g- t k = - ( 1 - XpdMf )e s . 

The pressure-dilatation covariance is more important than the dilatational dissipation and will act 
to slow the rate of decrease of k by shunting energy from the internal modes (mean temperature) 
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where it is stored into the kinetic modes. The results of the present analysis and Sarkar et al. 
(1991) are both internally consistent: they treat two different problems. _ The present analysis 
treats a turbulence in which the compressible portions of the flow are generated by the turbulent 
motions. Sarkar et al (1991) treat a turbulence on which is superposed, by the initial conditions, an 
Mf compressible velocity field. The postulated M? initial condition gives rise to an Mf dilatational 
field. The effects of compressibility will reflect the evolution of the compressible field in the presence 
of an incompressible turbulence and the dilatational field will then be order causing an increase 
in the k decay. The dilatational field is an order Mf effect when the dilatational fluctuations are 
generated, not by the initial conditions, but by the vortical fluctuations. These very interesting and 
potentially contradictory issues need to be investigated more closely with a DNS in which special 
care is taken with the implementation of the initial conditions. 

o 

Consider now a near equilibrium, T ~ 0, homogeneous, T* = Tj = 0, high Reynolds number, 
Pt — 0 shear flow. In such a flow the pressure-dilatation is now 

3 

< pd >= - XpdM ? [P k - <p>e- -M? 7(7 - 1) < p > e] . (97) 

There are several things worth noting. The first is the change of sign of < pd > noted by Sarkar et 
al. (1991a). For flows with small turbulence production the pressure dilatation is positive. If the 
production is large and exceeds dissipation by a certain amount the pressure-dilatation covariance 
is negative and there is a net transfer of energy from the turbulence field to the mean internal 
energy. This is consistent with numerical results of Sarkar et al. (1991a) only here the analysis 
indicates when the change of sign of < pd > occurs. 

O O 

The T and k equations for the homogeneous shear flow can be written 

§i k = (1 - X pdM?)[P k -e] + ^XpdMf 7(7 ~l)f (98) 

c vj^ T = XpdM? [P k - e] + e - |x P d-^< 4 7(7 - l)e (99) 

and the dependence on the production-dissipation balance and the the turbulent Mach number 
are more readily seen. Note that XpdM? < 1. Thus, to lowest order, 0(M?), the effects of the 
compressible dilatation is to reduce the excess production over dissipation in the k equation by a 
factor XpdM ? ; this energy is transferred to internal modes and the rate of increase of the mean 
temperature is amplified by an additive factor XpdM? [P k - e]. The reduction of the turbulence 
related quantities seen in compressible flows appears to be attributable to the transfer of energy 
from kinetic to internal modes. Computations with this pressure-dilatation representation in the 
compressible mixing layer have shown it to be primarily responsible for the substantial reduction 
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in growth rates. Additional contributions come from the T relaxations! term not carried in the 
equations above. These computations, as they represent the use of these representations for practi- 
cal calculations, are the subject of articles now in progress whose objective is to evolve the present 
mathematical results into a computational model. 


If in the representation for < pd > one sets < pd >= 0 one can predict the critical M ? c , for a given 
shear rate, turbulent Reynolds number and anisotropy, at which the pressure-dilatation changes 
sign. Below this critical turbulent Mach number there is a net transfer of energy from the mean 
temperature field in which it is stored. Above this critical turbulent Mach number there is a transfer 
and storage of energy in the temperature of the fluid rather than increasing the kinetic energy of 
the turbulence. Using Pk = \b12kS and Pj = 2 < p > S 2 the critical Mach number as a function 
of the anisotropy and shear rate can be obtained. To zeroeth-order 


2 _ 4 \ k 612 S- < p > ( 

1 3 7(7 - 1)[2 < p > S 2 + < p > (] 


( 100 ) 


Recall that e = e s + € c and that in fact the equation is a quadratic for A/ 2 if the dilatational 
dissipations contribution is included. If the dilatational dissipation is not distinguished a little 
more algebraic manipulation produces 


M{ = 


\ bn S - 1 


( 101 ) 


3 7 (7-l)[l + ^5 2 ] 

This is essentially a statement of the fact that as long as the production exceeds the dissipation by 
and order M ? quantity 


Pk > [1 + -M t 2 7(7 - 1)] <p>e 


( 102 ) 


that there is a net transfer and storage of energy in the mean temperature field. Which immediately 
suggests a numerical simulation of a homogeneous shear with an isotropic initial condition. At a 
critical Mach number < pd > will change sign as a function of the anisotropy and Reynolds number 
of the flow. 


5. Discussion and clarification of limitations and assumptions 

A few assumptions have been made to obtain representations for the dilatational covariances. The 
assumptions are not in anyway unreasonable but do limit the applicability of the results to specific 
classes of flows. This section is a compendium of the assumptions; it exists in order to insure that 
the applications of these representations be made with an awareness of their limitations. It should 
also help asses how much of the physics these representations capture or neglect in any specific 
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flow. It is hoped that full disclosure of the assumptions will suggest future work to account for 
their potential shortcomings. 

1) To model the effects of compressibility on the turbulence it has been assumed that the turbulent 
Mach number of the fluctuations, Mf — is small. This appears to be the case for several classes 
of supersonic flows of current engineering interest. The expressions derived are not expected to be 
useful for hypersonic flows in which eddy shocklets are important; such a flow situation will, in all 
likelihood, require a very different analysis. 

As has been discussed, the low M t assumption is equivalent to the compact source assumption 
of aeroacoustics. In the present context, this means that the correlation length scale of the flow 
structures producing the dilatational field is much smaller than the wavelength of the propagating 
field produced. This allows the dilatation to be algebraically related to the instantaneous material 
derivatives of the pressure fluctuations of the solenoidal field. 

The low turbulent Mach number assumption should not be understood to imply a low mean flow 
Mach number. 

It should also not, necessarily, imply a low gradient Mach number. The leading order contribu- 
tion to the Reynolds stresses, in the low r M t limit, is from the solenoidal field and the lack of 
signal communication across an eddy, when the gradient Mach number is high, effects higher order 
corrections to the Reynolds stresses. The present analysis treats the leading order dilatational 
fluctuations. The form of the equations derived suggests a ’thermodynamic 1 rather than signal 
propagation interpretation of the the gradient Mach number; the gradient Mach number scales 
the relative magnitude of the dilatational fluctuations to the pressure fluctuations generated by the 
mean velocity gradients. Interpreting it as a quantity characterizing the propagation of information 
is expected to be important in large turbulent Mach number situations. 

2) The scalings employed imply that the equations derived do not account for effects associated 
with phenomena that have coherences on much larger length scales. This would include variances 
and covariances of quantities that propagate. Such is the case with dilatation which has propagated 
into the local turbulence volume from regions more distant than the local integral scale. Though 
the correlation with the local turbulence is expected to be negligible, such signals are correlated 
with themselves and will make a contribution to the local dilatational dissipation. These effects 
are not accounted for in the present development; as such the development is limited to turbulence 
fields that are on order of or smaller than the acoustic scale of the flow. This is the compact flow 
assumption, D/\ < 1, and constitutes a limitation to the current representation. The limitation 
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may be nominal for many flows of current engineering interest in which the compact flow assumption 
is expected to be useful. 

The compact flow assumption also excludes the cumulative effects, over scales large with respect 
to the wavelength, the propagating field on the turbulence, as sometimes occurs in the sound 
generation and propagation problem. 

3) The covariances of the fluctuating dilatational field are assumed to be generated by and evolve in 
accordance with the turbulence field. They do not result from any externally imposed “acoustical” 
fields, or radiation from far field turbulence. Nor are they an adjustment of the flow to initial 
conditions with an acoustical or compressible component not generated by, or otherwise unrelated 
to the turbulent flow. 

The scalings and analysis employed imply that the source of compressibility in the flow is due to 
the turbulence within an integral scale of the position in question. The effects of q 2 /q 2 ~ 0(M f) 
compressible velocity fields, as seen in Sarkar et al (1991), superposed on the flow T by the initial 
conditions will produce much higher dilatational dissipation rates. This sensitive dependence of the 
dilatatio.nal dissipation rates on the initial conditions has been seen in the simulations of Blaisdell 
et al (1993). There are, in all likelihood, complex flow r situations in which sizable compressible 
fields, q 2 /q 2 are generated. This, however, is not the situation for which the present pseudo-sound 
theory has been developed. 

Here compressibility effects due to the finite Mach number of the vortical fluctuations are studied. 
IN studies of such compressibility effects resulting from the turbulence appropriate initial conditions 
are required for DNS. Consistency requires that the density and temperature variances be related 
to the incompressible field according to 7 < pp >= < 96 > = < pp >. Where, in dimensional 

terms, for an isotropic turbulence with no mean velocity gradients the pressure variance is 

<PP> S = \plok 2 Il (103) 

w'hile for a turbulence with divergence free mean velocity gradients, to leading order, 

< PP > T = ^Pl jt 2 [3S 2 + 5 W 2 } I[ = ^PL(|) 4 ^ 2 [35 2 + 5 W 2 } /[. (104) 

Or, in the context of a DNS starting from incompressible initial conditions, a point- wise propor- 
tionality between the fluctuating pressure and the density and temperature is required. The proper 
initial condition on the dilatational field is more difficult but may be much less important. 

4) The mean pressure and mean density has been assumed locally constant - constant over a length 
scale over W'hich the turbulence is correlated. This is also equivalent to the statement that the 
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sound speed is also locally constant; the stochastic nature of the fluctuations in the sound speed 
has been neglected. A corollary to the locally constant mean density assumption, is the fact that 
the mean dilatation is negligible. Thus, for flows in which the mean density and mean pressure 
vary appreciably over an integral scale the present representation is only a zeroeth-order theory. 

5) The turbulent fluctuations that contribute to the dilatational covariances have been assumed to 
obey the adiabatic gas law. This uncouples the problem from noisentropic aspects of a compressible 
flow which may be important in wall bounded flows. A scale analysis of the fluctuating dilatation in 
the near w r all region does appear to indicate that the nonisentropic contributions are higher order. 
This may not be the case for walls with large heat transfer. 

6) All the expressions presented have been obtained assuming that the major contribution to the 
quantities come from the isotropic portions of the statistics of the fluctuating field. The expres- 
sions obtained are the lowest order expressions in a series expansion in powers of the single-point 
anisotropy of the turbulence. Higher order terms allowing for contributions from the anisotropy 
are straightforward in concept but complicated in execution. This has been found to be the case 
in a few cases for which such parameterizations have been worked out. Perhaps, if the results are 
to be used for engineering calculations, the constants derived should be viewed as requiring some 
modifications for the anisotropy of the flow. A DNS might be useful to see which of the several 
contributions to the dilatational covariances are most sensitive on the anisotropy of the turbulence 
field. 

7) Throughout the development homogeneity has been invoked to make the statistical manipulation 
tractable. In practice this requires that the turbulent field be homogeneous on a scale 1/ L < 1 
where L is the scale of the inhomogeneity. Clearly few engineering flows meet this requirement; 
however, any representation that is created must at least be consistent with results obtained using 
this state. 

8) The assumption of quasi-normality has been made to achieve the statistical closure for the large 
scales of the flow. This involves the neglect of third-order moments with respect to second and 
fourth moments. This is an assumption that is extensively used and discussed throughout the 
literature. Corrections to the derived relations including the third-order moments are thought to 
be minor. 

9) The spectrum of the turbulence is assumed to have a negative power law behavior implying 
that a spectral Mach number of the fluctuations decays with wave number. This suggests that 
the portion of the spectrum exhibiting compressibility effects are at the lower wave numbers. It 
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is for this reason that the viscous effects have not been carried in the analysis for the dilatational 
covariances. Such effects may well need to be incorporated in near wall flows. 

10) The scale relation e = a(^) 3 / 2 /f has been used several times. While having substantial 
empirical support in incompressible flows it is now being used in situations in which its validity 
must be assessed. Is 6 = a(y) 3 / 2 /^ valid for a compressible turbulence ? If not how must it be 
modified ? What is the best definition of a length scale with which to define a gradient Mach 
number ? Implicit in the present derivations has been a longitudinal length scale. Sarkar (1994) 
uses a transverse length scale; how does this compare to the longitudinal length scale of an isotropic 
turbulence. Should c in the length scale definition include the compressible dilatation e c ? Clearly 
since the model is going to be used in sheared flows and near walls one must asses how accurate 
this relationship is for such situations. Sreenivasan (1984, 1994) has addressed the effect of shear 
on the relationship. 

6. Summary and Conclusions 

The mathematical consequences of a few assumptions about the size of the fluctuating pressure 
and density in a compressible turbulence are followed. A low turbulent Mach number singular per- 
turbation has produced a diagnostic constitutive relationship relating the fluctuating dilatation to 
the fluctuating pressure and velocity fields. This constitutive relation is the lynchpin of the devel- 
opment allowing closure for the effects of compressibility in terms of the divergence-free portions of 
the fluctuating flow field. Moments of the constitutive relation produce analytically consistent rep- 
resentations for the dilatation variance and the pressure-dilatation covariance in a turbulence field 
with and without mean velocity gradients. Application of the methods of statistical fluid mechanics 
and the assumptions of quasi- homogeneity, quasi-normal behavior, and isotropy produces expres- 
sions for the covariances with the fluctuating dilatation. Except for the well-established empirical 
result, l ~ (|fc) 3 ^ 2 /c 5 , used to close the expressions, no additional phenomenological assumptions 
are made. The analysis is, in the low limit, exact and produces representations for the effects of 
compressibility in which there are no undefined constants. The constants that appear are known in 
terms of integrals of the longitudinal velocity correlation of an incompressible isotropic turbulence. 

Both Lele (1994) and Blaisdell et al (1994) have reflected that an algebraic closure for the effects 
of compressibility solely dependent on M 2 appears to be overly restrictive. The present analysis 
has quite naturally indicated the importance several additional parameters. The compressible 
dissipation is found to be a function of the local values of the turbulent Reynolds number, R t * 
the turbulent Mach number, Af#, the two mean velocity gradient Mach numbers Ms and M\y - 
and the solenoidal dissipation. The pressure-dilatation is seen to be a nonequilibrium phenomena. 
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It is found to be a function of the rate of change of the turbulent Mach number, M t , the mean 
density, < p >, the energy of the turbulence, k , and the two relative time scales Sk/e s and 
Wk/e s . Additional manipulations show that it can be expressed as a function of the production, 
the dissipation and the transport. There still remain the nonequilibrium effects associated with the 
adjustment time of the eddy turnover time to the time scale of the mean flow. 

The analysis has produced a few simple and interesting metaphors for the effects of compressibility 
on a turbulent flows. For classes of turbulent flows for which this analysis is relevant, one example 
being the compressible mixing layer, the results suggest mechanisms that play a role in the reduction 
of the mixing layer growth rate. In short, the pressure-dilatation transfers turbulent kinetic energy 
to the internal energy field effectively reducing the relative excess of production over dissipation by 
a factor XpdM Further reduction of k occurs through the dilatational dissipation. These effects 
are dependent on the local mean flow gradients and as such suppress growth rates most in regions 
of high production. 

It is hoped that these results will be of use in further understanding the complex effects of com- 
pressibility and stimulate additional new investigations including the assessment of the consequences 
and utility of the present pseudo-sound development. Some of the results worth noting are now 
summarized: 

1) Noteworthy is the appearance of two mean flow Mach numbers based on the mean velocity 
gradients. One based on the mean deformation, Ms = Sl/c , and another on the mean rotation, 
Mw = Wife. These gradient Mach numbers have been identified as important parameters in 
assessing the effects of compressibility in the numerical experiments of Sarkar (1994). 

2) The pressure-dilatation is essentially a reversible nonequilibrium phenomenon acting as a mech- 
anism by which the fluctuating kinetic energy is transferred and stored in the internal energy field. 
The pressure-dilatation is shown to be interpretable as an added mass effect reducing the rate of 
change of turbulence quantities by the capacitance of the mean internal energy field. As such it 
may, in part, be responsible for some of the relaxational effects seen in compressible flows. This 
remains to be seen. 

3) The rapid portion of the pressure-dilatation is seen to be a function of the rate of change of 
the relative times scales, Sk/e s and Wk/e s . These quantities may be interpreted as indicating 
how closely the eddv-turnover time tracks the mean velocity gradient. They occur because of the 
evolution of the length scale of the turbulence and its importance to the rapid pressure integral. 
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4) Unlike the solenoidal dissipation, the theory predicts that the compressible dilatation is a function 
of the Reynolds number and vanishes, for fixed turbulent Mach number, as. the Reynolds number 
increase. The compressible dilatation, in the low M t limit, doesn’t appear to be interpretable as a 
spectral flux. 

A brief overview of the more specific results of the theory can be found in the introductory section. 
Limitations and extensions of the theory have been indicated in the previous section. A summary 
of the representations is given in the appendix. 

This article is a documentation of the physical implications of a pseudo-sound analysis for for the 
covariances of the fluctuating dilatation. The uniqueness of the investigation is the small number 
of phenomenological assumptions made. The results are a mathematical consequence of the initial 
assumptions. As an analytical investigation the article is complete. The assessment and utility of 
these results as models for engineering computations is the subject of several works planned or in 
progress, Ristorcelli (1995), Ristorcelli et al (1995). Preliminary computations in a few r simple 
benchmark flows have been successful. 
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Appendix 1: Synopsis of the dilatational covariance representations 

The analytical results, and associated mathematical formulae are, as an aid to the user, briefly 
summarized. It should be emphasized that, at this point, the representations do not constitute a 
fully developed or tested turbulence model. Testing, verifying and evolving the present analytical 
results into a working turbulence model suitable for engineering calculations is the subject of work 
now in progress. 

The kinetic energy equation for the turbulence is 

<p>j^k = Pk-<p>€ s + <pd>-<p>€ c + Tk ( 105) 

where Pk represents the production and T* represents the transport terms. The effects of compress- 
ibility are given by the second two terms. The dissipation e s = u < c OjU 3 > is the usual dissipation 
associated with the vortical motions of the incompressible turbulence. The dilatational dissipation 
or the compressible dissipation and denoted e c = < dd >. There are additional terms repre- 

senting the the contraction of the dyad of the mass flux on the mean flow acceleration. T will be 
used to represent all such terms as well. Local isotropy has been assumed for the dissipation. The 
usual modeled dissipation equation 

D/Dt e = - ( C £ i < UiUj > Ui,j + C e2 £) e/k + T t (106) 

is carried to describe its evolution. Note that no corrections for compressibility have been made in 
this equation. The mean temperature equations , with the assumption of constant c v , is 

< p > Cyjft T = Pj — <pd> + <p> e 5 + < P > e c + Tj 

Where Tj is the transport of the mean temperature including such effects as the mean heat flux 
and the turbulent or pressure transport. The production for a homogeneous flow with homogeneous 
mean velocity gradients is Pj = -PD + 2 < p > (77i, 2 ) 2 - Here P is the mean pressure and D 
is the mean dilatation. Depending on the particular application the temperature equation may 
have more or fewer production terms than indicated here. Note that if one carries the total energy 
equation in a simulation, rather than the mean temperature, that that the pressure-dilatation only 
needs to be carried in the k equation. 


The representations for the effects of the compressible dissipation are given by the sum of the slow 
and rapid portions, e c = e r c + ef. 
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M t is the turbulent Mach number, M 2 = | k/c 2 , where = 7 < p > / < p > is the local sound 
speed. Note that its definition follows that used in the acoustic literature where such a Mach number 
is traditionally used to describe the sound generation problem by turbulence. It is the parameter 
that arises naturally in nondimensionalization of the equations for a small parameter expansion. 
The turbulent Reynolds number is given by R t = ^ = using the facts that u = 2Ar/3 and e s = 
au 3 /£ which is used to express the length scale as, l = a(2k/3) 3 ^ 2 /e s . For simple incompressible 
shear flows, the constant, a, varies between 1 - 4 depending on flow, Sreenivasan (1995). Note that 
in the definition the characteristic velocity | k is used; not < U\U\ > as is sometimes the case. The 
nondimensional strain and rotation rates are given by: S 2 = ( Sk/e s ) 2 , W 2 = (Wk/€ s ) 2 where of 
course, 5 = \/Si 3 S t 3 and W — \JW X3 Wi 3 . The strain and rotation tensors are defined in analogy 
with the incompressible case. ie. traceless S ZJ = +Uj*i — W ZJ = ~U JU ]. Note 

that Sjj = 0 since D — f/j, r In a simple planar shear flow, ?7;,j = U 1,2 that S 2 — W 2 = |f r2 ,2- 
A quick of order of magnitude estimate for the integrals can be made using / = e“^ 2?r / 4 . The 
following values are found: If = ^, /| = — 4.77, JJ = f = 0.785, /[ = £ = 1.273, = 3. 

The values found from high Reynolds number wind tunnel data are different: If = 0.300, /f = 
13.768, /j = 2.623, If = 1.392, = 3, Zhou (1995). The values given for the integrals reflect the 

assumption of an equilibrium isotropic turbulence and are to be understood as suggestive of the 
order of magnitude that they may have in more complex anisotropic and inhomogeneous situations. 

The full pressure-dilatation covariance is a sum of two terms, < pd > — < pd > r 4- < pd > 5 , 

< Pd > s = - §- t i<P> M?k] 

< pd > r = - ^ (|) 3 /f a 2 ~ [< p > kM 2 [3S 2 + 5 W 2 } ] (108) 

Summing and using the evolution equations for T and k produces the following quasi-algebraic 
representation for the full pressure-dilatation covariance, 

< pd >= -x P dM? [ Pk - < p > € + T k - f M, 2 7(7 - 1 )(Pt+ < P > e + T t )\ 

- < p > k M 2 x r pd &T 

2 Lpd 

1 + 2I p dM 2 + - 1) 

& 

i + 2i pd M? + - i ) 

\n + iu [ 3^ 2 + bw 2 ] 

^(| )3 Q 2 A r - (109) 


Xpd ~ 

X r pd = 
lpd — 
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Note that e = e s + e c and T = [35 2 + 5 W 2 ]. The term inside the inner brackets is the right hand 
side of the mean temperature equation. 
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